Ideal strength of ferromagnetic Fe-based alloys from first-principles theory 
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The all-electron exact mufRn-tin orbitals method in combination with the coherent-potential ap- 
proximation has been employed to investigate the ideal tensile strengths of bcc ferromagnetic Fe and 
Fei-xMx (M=Cr, Ni, Al, Co, Mn, and V) random alloys in the [001] direction. The present ideal 
strength is calculated to be 11.0 GPa for pure bcc Fe, in good agreement with the available theoret- 
ical data. For the Fe-based alloys, we predict that Co, Cr and V increase the ideal tensile strength, 
while Ni and Al decrease it. Manganese yields a non- monotonous alloying behavior. We show that 
the limited use of the previously established ideal tensile strengths model based on structural energy 
differences in the case of Fe-bases alloys is attributed to the effect of magnetism. 
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I. INTRODUCTION 

The mechanical strength of real materials is usually 
controlled by the occurrence of grain boundaries, cracks, 
dislocations and other micro-structural defects. If such 
defects were not present, the strength would be limited 
by the stress at which the lattice itself becomes unsta- 
ble with respect to a homogeneous strain. This stress 
referred to as the ideal strength or ideal tensile strength 
(ITS) is an upper bound strength of an ideal single crys- 
tal. The ideal strength is an inherent property of a defect- 
free material and can offer insight into the correlation 
between the intrinsic chemical bonding and the crystal 
symmetry, and has been accepted as an essential me- 
chanical parameter of single crystal materials."' The ex- 
perimental data on the ideal strength are rather limited. 
The few available values were obtained from tensile tests 
for whiskers^'^ and from nanoindentation experiments.** 

In recent years, considerable attention was paid to the 
computation of the ITS of elements or ordered alloys by 
first-principle methods. The ideal strength of refractory 
metals (such as Mo, Nb, V, and W), noble transition 
metals (such as Cu, Pt and Au), elements crystallizing 
in the diamond structure (Si, Ge, and C) as well as few 
ordered alloys (TiAl and NiaAl) was extensively investi- 
gated.^"-'^ Among the body-centered cubic (bcc) metals, 
Fe received a great scientific interest. '''^"^' It was re- 
ported that, the [001] direction of Fe is the weakest direc- 
tion in response to an uniaxial stress. It has been shown, 
that a bifurcation from the primary tetragonal deforma- 
tion path (the Bain path) to a secondary orthorhombic 
deformation path occurs after the ITS on the primary 
path is reached.^ Thus, bcc Fe strained along the [001] 
direction fails by cleavage and not by shear as opposed to 



bcc V and Nb.''^ In spite of all these theoretical efforts, 
the first-principles description of the ITS in substitution- 
ally disordered alloys is rather limited. Li et al.^^ inves- 
tigated the composition dependent ITS of bcc vanadium- 
based nonmagnetic (NM) random solid solutions alloyed 
with chromium and titanium in various crystallographic 
directions employing the coherent potential approxima- 
tion (CPA). To the best of our knowledge, no experi- 
mental or theoretical study has focused on the alloying 
effect on the ITS of Fe-based alloys. It is the purpose 
of this paper to bring forward a comprehensive study of 
the ideal strength to bcc ferromagnetic (FM) Fe-based 
random alloys. 

Iron alloys are the most widely used engineering ma- 
terials due to their excellent mechanical properties. For 
example ferritic Fe-Cr-based stainless steels have been 
considered as the primary structural materials in the first 
wall and blanket structure of future fusion reactor. ^■^■^'^ 
Alloying plays a central role in designing advanced engi- 
neering materials with desired properties. Different so- 
lute atoms produce different effects on fundamental me- 
chanical properties. Previous works focused on the ef- 
fects of various typical solute atoms on the mechanical 
properties of Fe in the small deformation region, where 
the stress-strain relations are linear. ^^"•^° However, the 
ideal strength describes the mechanical properties of the 
material beyond the elastic region. 

In this work, using the exact muffin-tin orbitals 
method, we study the composition dependence of the 
ITS of FM ¥ei-xMx alloys for various alloying elements, 
M=Cr, Ni, Al, Co, Mn, and V. The primary purpose 
of our work is to give an account of the application of 
the alloy theory to the ITS of bcc FM random alloys. 
Second, we aim to provide a consistent theoretical guide 



to further optimization of the composition of Fe-based 
alloys in multiscale materials design. 

The structure of the manuscript is as follows. In Sec- 
tion II, we describe the computation tool and the im- 
portant numerical details. The results are presented and 
discussed in Section III. Here, first we assess the accu- 
racy of our calculations by considering pure bcc Fe. Then 
we study the effects of the alloying elements on the ideal 
strength of Fe-based alloys. 



II. COMPUTATIONAL METHOD 
A. Total energy calculation 

The first-principle method used in this work is based on 
density functional theory (DFT).^"'^ We adopted the gen- 
eralized gradient approximation (GGA) of the Perdew- 
Burke-Ernzerhof (PBE)"^^ functional to describe the 
exchange-correlation interaction, which is well known 
to give the correct FM bcc ground state for Fe. The 
Kohn-Sham equations were solved using the EMTO 
method. ■''^^■'''' The problem of disorder was treated within 
the CPA and the total energy is computed via the full 
charge-density technique. '^^"^^ 

The EMTO method is an improved screened Korringa- 
Kohn-Rostoker (KKR) method, ^'^ where the full poten- 
tial is represented by overlapping muffin-tin potential 
spheres. Inside these spheres, the potential is spheri- 
cally symmetric and constant in between. By using over- 
lapping spheres, one describes more accurately the ex- 
act crystal potential compared to conventional muffin-tin 
methods. Further details about the EMTO method and 
its self-consistent implementation can be found in previ- 
ous works. ^^"^^'^^'^^ The accuracy of the EMTO method 
for the equation of state, elastic properties, and the ideal 
strength in tension of metals and alloys has been demon- 
strated in a number of previous works. ^^■'^^'^^'''""^^ 

The paramagnetic state of the Fe-based alloys was sim- 
ulated by the so-called disordered local moment (DLM) 
model. ''^ Within the DLM picture, the paramagnetic Fe 
and Fei^xMx binary alloys were described as a binary 
alloy FetFel and a quaternary (FetFe|)i_j;(M fAf i)x 
alloy, with an equal amount of spin-up (^) and spin-down 
(t) components, respectively. 



B. Ideal tensile strength calculations for bcc 
crystals 

The principles of the response of bcc crystals to uni- 
axial loading were developed in a series of works by Mil- 
stein et al.*''^^*^ In this work, we computed the ITS in 
the [001] direction of bcc Fe and Fe-based alloys. Since 
[001] was already identified to be the weakest direction 
of bcc Fe,"'^*'^^'^^ here we concentrated on this direction 
only. Assuming an uniaxial tensile load, the tensile stress 



a{e) can be calculated by 



a(e) 



15 



1 dE 



(1) 



where E is the total energy per atom and r2(e) is the vol- 
ume per atom at a given tensile strain, e. The first max- 
imum on the stress-strain curve defines the ITS, am, for 
the selected deformation path with corresponding maxi- 
mum strain, em. 

Uniaxial loading along the [001] direction reduces the 
symmetry of the bcc lattice to the body-centered tetrag- 
onal (bet) one on the primary deformation path. Clat- 
terbuck et al} reported a bifurcation from the primary 
deformation path to a secondary orthorhombic (orth) de- 
formation path in Fe, however, the branching occurs for 
strains well above em- Thus, an isotropic Poisson con- 
traction along the Bain transformation path describes 
appropriately the symmetry of the distorted Fe lattice 
up to £„!. The ideal tensile strength of the present Fe- 
alloys may nevertheless be influenced if this branching 
point shifts towards strains smaller than the maximum 
strain along the primary deformation path. We account 
for this possibility in this work. 

On the primary tetragonal deformation path, the uni- 
axial strain energy,^^ AE{c\ [001]), describes the total en- 
ergy change upon deforming the bcc crystal in the [001] 
direction, and relaxing it with respect to the dimensions 
in the (001) plane. 



A.B(c; [001]) = min E{a, c) - Eq 



(2) 



here a and c denote the lattice constants of the quadratic 
basal plane and the height of the bet unit cell, respec- 
tively. The initial undistorted state corresponds to the 
equilibrium bcc structure with energy Eq and c — a. At 
c/a — \/2, the bet lattice coincides with the face-centered 
cubic (fee) lattice. On the secondary orthorhombic de- 
formation path, we consider A£^(aorth, ^orth; [001]) and 
minimize the total energy with respect to the lattice pa- 
rameters of the face-centered orthorhombic lattice, Corth 
and 6orth (to describe the bifurcation from the primary 
tetragonal to the secondary orthorhombic strain path, 
the face centered tetragonal reference frame of the Bain 
transformation is used^^). 

For a more detailed technical description of ideal 
strength simulations with EMTO, we also refer the reader 
to Ref. 22 



III. RESULTS 

A. Ideal strength of bcc iron 

To assess the reliability of our computational approach, 
we first performed the simulation of a tensile test in 
FM bcc Fe for uniaxial loading along the [001] direction. 
The theoretical equilibrium lattice parameter of bcc Fe 



is 2.829 A from our calculation, which agrees well with a 
full-potential value of, 2.830 A.^° The experimental lat- 
tice parameter, 2.868 A,^^ is 1.4% larger than the theo- 
retical equilibrium value. The calculated spin magnetic 
moment of FM bcc Fc is 2.27 fiB, which is approximately 
2 % larger than the measured value 2.22 fiB-^^ This over- 
estimation of the magnetic moment of Fe is due to the 
employed GGA functional. 

The magnetic order of the strained Fe lattice may 
change along the uniaxial loading process, which is for 
strains up to the branching point governed by the Bain 
transformation as discussed in the previous section. At 
zero strain, the ground state magnetic order is FM (bcc 
phase) . The fee state of Fe lying on the Bain transforma- 
tion path (albeit at strains much larger than em) exhibits 
a non-coUinear spin arrangement as measured in fee Fe 
precipitates and in thin fee Fe films. ^^"^^ That indicates 
that the ground state magnetic order of the strained bet 
Fe lattice begins to differ from FM order at some partic- 
ular strain in the range between bcc and fee along the 
primary transformation path. If this would be the case 
for e < em, then additional magnetic order should be 
considered. However, there are strong indications that 
the FM order is the prevailing magnetic state for strains 
smaller than and somewhat above em.^°'^^'^^ 

Clatterbuck et al.^'^ computed the ideal strength of 
Fe in the [001] direction considering the FM order, a 
coUinear anti-ferromagnetic structure (AFM, magnetic 
moment sequence fi) and a coUinear double layer anti- 
ferromagnetic structure (DAFM, magnetic moment se- 
quence ttii)- Their results show that Fe remains FM 
up to the point of its elastic instability during uniax- 
ial tension, which lies at approximately em — 15 % with 
c/a « 1.16 and fl/flcxp > 1 (^oxp denotes the experi- 
mental atomic volume of bcc Fe). Further evidence is 
given by Tsetseris^® and Friak et al. which published 
minimum energy contour plots with respect to various 
magnetic orders (FM, AFM, and DAFM order by Friak 
et al. and non-coUinear magnetism via a spin spiral for- 
malism by Tsetseris) as a function of the bet geometry 
thereby defining magnetic phase boundaries between dif- 
ferent magnetic states. According to both references, FM 
order is the predominant magnetic order in the configu- 
ration space for c/a < 1.25 and il/ilexp > 0.95. Both ref- 
erences hence indicate that the point of elastic instability 
of Fe reported by Clatterbuck et al. is indeed located far 
from the borderline of FM order towards any other in- 
vestigated magnetic ground state order. Hence, here we 
assume that Fe and also the present Fc-rich binary alloys 
stay FM during the deformation process. 

The ideal tensile strength a^ corresponding to the 
strain em from our and other calculations in [001] di- 
rection for iron are listed in Table I. It can be seen that 
the three projector-augmented wave (PAW) works^'^^'^^ 
show similar stress values, however scatter somewhat in 
em- Namely, Clatter et al.^ reported a value of a^ = 12.6 
GPa at £m = 15% and Liu et al}^ gave a value of 
(Tm=12.4 GPa at em=14%. Furthermore, the two cTm 



TABLE I. The ideal tensile strength (Jm and the correspond- 
ing strain Em from our and published computations in [001] 
direction for a pure iron crystal. FP-LAPW stands for full- 
potential linearized augmented plane wave. 





method 




direction/ [001] 




an. (GPa) 


em{%) 




EMTO 




11.0 


14 




PAW Ref. 1 




12.6 


15 


Fe 


PAW Ref. 18 




12.4 


16 


PAW Ref. 19 




12.4 


14 




FP-LAPW Ref. 


20 


14.2 


15 




FP-LAPW Ref. 


21 


12.7 


15 



values obtained by the all-electron full-potential linear 
augmented plane wave (FP-LAPW) method differ by 
1.5 GPa {^ 11%). Compared to these theoretical full- 
potential data, the EMTO result for the stress of Fe is 
slightly smaller but still in reasonably good agreement. 
We obtained c/a « 1.18 and fl/ft 



cxp 



1.04 for the 

point of instability of Fe in agreement with Clatterbuck 
et o/.^" Finally, we computed that the bet to orth branch- 
ing occurs at e = 19 % well above em in accordance with 
Ref. 1. 

On the experimental side, the only attempt to mea- 
sure the ideal tensile strength of Fe was by Brenner^ 
who tested Fe whiskers for tension was measured to be 
the [001] direction. The ITS of these Fe whiskers at 
room temperature was measured to be approximately 5 
GPa. Compared to the theoretical data, the experimen- 
tal value seems to be considerably lower. However, the 
observed failure initiated at the surface and, therefore, 
the measured value cannot be considered to represent 
bulk strength. 

Based on the above results, we concluded that our the- 
oretical tool is able to describe the ITS of bcc Fe with 
sufficiently high accuracy. 



B. Ideal strength of Fe- based alloys 

In the following we turn to bee Fe-based alloys and in- 
vestigate the effect of alloying elements M=Cr, Ni, Al, 
Co, Mn, and V on the ITS of Fei_a;Ma; random solid so- 
lutions. The selected solute atoms are common in com- 
mercial Fe-based steel alloys and they represent simple 
metal (Al), nonmagnetic (V) and magnetic (Cr, Co, Ni, 
and Mn) transition metals. The total concentration of 
these solutes was varied in the range from to 10% 
except for Mn which the maximum concentration was 
5%. The theoretical equilibrium lattice parameters of 
Fei-xM^ (M=Cr, Ni, Al, Co, Mn and V) random alloys 
are displayed in Table II. We can see that compared to 
pure Fe, all alloying elements increase the lattice con- 
stant. Our computed lattice constants practically repro- 
duce the earlier results from Ref. 28. 



TABLE II. Theoretical (EMTO) lattice parameters aix) (in A) calculated for the ferromagnetic bcc Fei_^Ma; alloys as a 
function of composition. For pure Fe, the value is a(0) = 2.829 A. 











a{x) 










Fe-Cr 


Fc-Co 


Fe-V 




Fe-Ni 


Fe-Al 


Fc-Mn 


2.5 


2.843 


2.837 


2.842 




2.843 


2.844 


2.839 


5 


2.849 


2.842 


2.849 




2.852 


2.853 


2.843 


7.5 


2.850 


2.845 


2.850 




2.858 


2.860 




10 


2.849 


2.848 


2.858 




2.861 


2.866 





Figure 1 shows the composition dependence of the ITS 
of binary bcc Fe-based random alloys along the [001] di- 
rection and the corresponding numerical data for selected 
compositions are listed in Table II. The calculated ideal 
strength is found to increase with Cr, Co and V and de- 
crease with Ni and Al addition to Fe. For instance, when 
10% Cr, Co, or V is added to bcc Fe, the ITS of Fe in- 
creases by 1.4, 0.8 and 2.4 GPa, respectively. If however 
10% Ni or Al are added to the Fe matrix, the ITS reduces 
by 1.5 GPa. For Mn, first the ITS increases weakly with 
increasing Mn amount below 2.5%, i.e., from 11.0 GPa to 
11.3 GPa, then reduces by 0.3 GPa when up to 5% Mn is 
added to Fe compared to pure Fe. According to Fig. 1, 
it is clear that Al, Ni, and V show the strongest average 
alloying effect (AcTm/Aa;) on the ITS of bcc Fe, whereas 
Mn, Ni and Cr produce intermediate alloying effect. 
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FIG. 1. (Color online) The ideal strength of Fe-based alloys 
as a function of composition. 



We also considered the possibility of a branching from 
the primary bet deformation path to the secondary or- 
thorhombic deformation path for the present Fe-based 
alloys. Based on these additional calculations, we can 
exclude a bifurcation from the primary bet deformation 
before the ideal strength is reached, i.e., for all alloys con- 
sidered here the branching occurs at strains larger than 



IV. DISCUSSIONS 



Failure of the structural energy diiference 
model for ITS 



Understanding the alloying effects on the ITS of Fe- 
based alloys is highly desirable. To this end, we start 
from a previously established simple model based in 
structural energy differences. For binary and ternary 
vanadium based Vi-y-zCryTi^ random solid solutions, 
0<J/ + 2<0.1, we established a qualitative correlation 
between the change of the ITS due to alloying and the 
change of the fcc-bcc structural energy difference (SED) 
of the alloy. ^^ Accordingly, the maximum stress {(J-i^ was 
approximated by cr^^ along the [001] direction defined 
as 



^SED 



1 A^SED 



Q. 



bcc 



(3) 



where AE^^^ denotes the SED at fixed volume (here 
of the bcc ground state, Jlbcc) and Ae is the strain at 
constant volume necessary to transform the bcc lattice 
into the fee lattice along the Bain transformation (Ae « 
0.260). 

Equation. (3) is based on the fact that the fee structure 
corresponds to the nearest maximum to the bcc phase of 
the uniaxial strain energy curve. ^'^^^^^ Thus, the uniaxial 
strain energy must level off to the fcc-bcc energy differ- 
ence, which implies a limitation on the maximum stress 
since it restricts AE{c; [001]) to the SED within the strain 
interval to accomplish the transformation from bcc to fee, 
i.e., the ratio AE^^^/Ae is bounded. 

The alloying trend is captured by Eq. (3) if there is 
a correlation between the change of the ITS as a func- 
tion of concentration, A<7m(a;), and the change of the 
fcc-bcc SED, A{AE^^^){x). The prefactor 1/rJbcc is 
weakly concentration dependent and does not change 
the conclusions drawn here. We show in the following 
that the correlation suggested by Eq. (3) fails for the 
bcc FM Fe-based alloys as opposed to the nonmagnetic 
Vi_j;^^Cr.yTi^ alloys.22 

First, we assumed a FM state for both fee Fe and the 
Fe-based alloys, since the FM fee state of Fe was shown 
to be the nearest maximum of the uniaxial strain energy 

„,,.^,,„ 1,19-21 



TABLE III. The ideal tensile strength am 


corresponding strain em from our calculations in [001] direction for iron-based alloys. 


o-m (GPa) 


em(%) am (GPa) £,„(%) a^ (GPa) £„,(%) 


"""'"'"""""" Fo-Cr 


Fe-Co Fc-V 



2.5 
5 

7.5 
10 

2.5 
5 

7.5 
10 



12.1 
12.4 
12.4 
12.4 

10.9 
10.6 
10.1 
9.5 



Fc-Ni 



14.0 
14.7 
14.5 
14.2 

14.3 
14.3 
13.6 
13.9 



11.1 
11.5 
11.7 
11.8 

10.5 
10.2 
9.7 
9.5 



Fe-Al 



14.2 
15.2 
14.9 
15.1 

14.7 
14.6 
14.4 
14.3 



11.9 
12.6 
13.2 

13.4 

11.3 
10.7 



Fe-Mn 



14.8 
15.3 
15.3 
15.3 

14.6 
14.4 



Figure 2 displays AfTm (x) , obtained from the values in 
Table III, as a function of A{AE^^^){x) calculated with 
EMTO-CPA. In figure, the alloying effect was obtained 
by increasing the concentration of the solute from % 
to 5 %. We can see that Ni, Mn, and Al decrease the 
ITS and also decrease the SED, however, Co, Cr, and V 
increase the ITS but decrease the SED. We also investi- 
gated the correlation between a^ and A£'^^^ increasing 
the concentration of the solute from 5 % to 10 %, however 
the result is qualitatively identical to the one depicted in 
Fig. 2. From these results, we infer that the anticipated 
correlation between the ITS and the FM SEDs following 
Eq. (3) can not capture the observed alloying trend of 
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FIG. 2. (Color online) The change of the ITS, Aam versus the 
change of the fcc-bcc SED, A{AE^^°), for Fei_^M^ alloys 
assuming a FM fee state for a concentration increase of the 
solute from 0% to 5%. Data in unshaded areas affirm a 
correlation based on Eq. (3). 



Next we tried to establish a correlation between 
the PM SEDs obtained from paramagnetic (PM) fee 
Fei-xMx alloys and their FM bee phases. Using the 
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FIG. 3. (Color online) The change of the ITS, Aam versus the 
change of the fcc-bcc SED, A{AE^^°), for Fei-a;M^ alloys 
assuming a paramagnetic fee state for a concentration increase 
of the solute from % to 5 %. Data in unshaded areas affirm 
a correlation based on Eq. (3). 



above described procedure and a concentration increase 
of the solute from % to 5 %, we plotted the stress change 
versus the PM SED change in Fig. 3. From that figure, it 
can be clearly seen that Ni and Mn decrease the ITS and 
likewise decrease the SED, while Cr, Co, and V increase 
the ITS as well as they increase the PM SED. However, 
FeAl is outlying as the ITS decreases with Al, but the 
PM SED increases with Al. Assuming a paramagnetic 
fee state thus captures somewhat better the observed al- 
loying trends than a FM fee but the correlation dictated 
by Eq. (3) is still not perfect. 

One possible reason why Eq. (3) is not a good estimate 
for the alloying trend of am for Fe alloys may be related 
to the different employed deformations paths. Along the 
unconstrained primary deformation path of tensile stress, 
the volume of the unit cell can change and a is deter- 
mined from Eq. (2). Fcc-bcc SEDs on the other hand 



were computed at the bcc equilibrium volume, i.e., on a 
Bain deformation path at constant volume ribco where a 
is given by 2ribcc/c^- Equation (3) captures the alloying 
trend on am, if the effect of alloying on the total energy 
is comparable for both strain paths. That was indeed the 
case for nonmagnetic V-based random alloys. ^^ However, 
for the present FM Fe-based alloys the interplay between 
magnetism and the structure in Fe along the deforma- 
tions paths should also be considered. We provide argu- 
ments in the following section, that the above failure can 
indeed be originated in the effect of magnetism. 



B. Magnetic effects on the ITS 

To achieve a better understanding of the matter, we 
computed an auxiliary ideal stress assuming a constant 
volume deformation, i.e., the maximum stress according 
to Eq. (1) was deduced from the total energy computed 
as a function of strain for fixed unit cell volume. The 
constant volume for each Fe-alloy is the one of its theo- 
retical equilibrium volume (Table II). The so-calculated 
maximum stress at constant volume, a^, with maximum 
strain, e^^, is then correlated to FM /\E^^^ from Eq. (3). 
Since both values were obtained for the same strain path, 
we expect that the correlation between ct^ and A£'^^^ 
significantly improves for all alloys, despite possible in- 
terference from magnetism. 

First, a^ as an approximate to the ideal strength of 
Fe-based alloys as a function of the alloying elements is 
displayed in Fig. 4. Notice the difference between Fig. 1 
and Fig. 4, the prior corresponding to relaxed volume and 
the latter to constrained volume. According to Fig. 4, the 
constrained ITS decreases with increasing composition 
for all alloying elements compared to the value of Fe. Al, 
Ni, and Mn solutes have a stronger effect than the other 
elements. 

Compared to the ideal stresses depicted in Fig. 1 (with 
variable volume), the constrained-ITSs keep their trends 
as a function of x in the case of Al, Ni, and Mn, but the 
ITS shows an opposite trend for Cr, Co, and V addition. 
Figure 5 displays the correlation between the change of 
the constrained-ITS calculated at the constant volume, 
Acr[,|(a;), and the change of the FM SED {/:i{/::^E^^^){x)). 
Both alloying effects for an increase of the concentration 
of the solute from % to 5 % and for an increase from 
5 % to 10 % arc displayed. Quite interestingly, we obtain 
that the change in the constrained-ITS correlates well 
with the SED change. 

By comparing Fig. 4 with Fig. 1, it is immediately ev- 
ident that the ideal strength of Fe is most significantly 
affected by the fixed volume constraint, i.e., it increases 
from 11.0 GPa to 16.0 GPa. This is unlike the Fe-alloys 
with moderate solute concentrations, where the ideal 
strength increases are considerably smaller the higher the 
solute concentration are, e.g., we obtained an increase 
from 13.4 GPa to 15.1 GPa for Fe-lOV. 

An answer to the observed correlation may be given 
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FIG. 4. (Color online) The ideal strength a^ of Fe-based 
alloys as a function of composition at a constant volume de- 
formation (closed symbols). The open pentagons are the ideal 
strengths of FeV alloys with fix magnetic moment along the 
relaxed volume deformation. 
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5. (Color online) The change of the ITS for constant vol- 
ume, AcrJ^ versus the change of the fcc-bcc SED, A(A_E^'^°), 
for Y^ei-xMx alloys assuming a FM fee state for a concen- 
tration increase of the solute from 5 % to 10 % and from % 
to 5 % (inset) . Data in unshaded areas affirm a correlation 
between Act," and /^E^^^ . 



by considering magnetism in Fe and its alloys. In the 
following we chose Fe-lOV as a representative of all in- 
vestigated binary alloys. The total energy and the spin 
magnetic moment of Fe are opposed to those of Fe-IOV 
in Fig. 6, where we plotted both quantities as a function 
of the tetragonal axial ratio {c/a) and the Wigner-Seitz 
radius, to. We also displayed the corresponding uniaxial 
deformation path in the range from e = . . . e,n. 

As visible from the contour plots, the magnetic mo- 
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FIG. 6. (Color online) Total energy (in mRy) and magnetic moments (in /ib) contours for pure Fe (upper panel) and V-lOV 
alloy (lower panel) as a function of the tetragonal ratio (c/a) and Wigner-Seitz radius (oj in bohrs). The energy are plotted 
relative to the bcc (cja — 1) minimum. Figures (a) and (c) are total energy contours, (b) and (d) are magnetic moments 
contours. The black solid lines show the uniaxial deformation paths in the range from e = ... em- 



inent increase of Fe along the deformation allowing for a 
volume change from e = . . . e,n is much stronger than 
the one of Fe-lOV; the numbers read 0.27 /ub and 0.09 ^b 
for Fe and the Fe-lOV alloy, respectively. These numbers 
are in contrast to the deformation at constant volume 
(in the range e — 0. . .e^J, where the changes are 0.09 
fiB ior Fe and 0.03 /^b for Fe-lOV. That indicates that 
the increase of the atomic volume has a more pronounced 
influence on the relaxation of the magnetic moment than 
the structural change {c/a), a result which can be clearly 
associated with the contour lines in Fig. 6. It should be 
noted that |e^ — e 
same material), i.e., e 
the present Fe-alloys. 

Now, alloying Fe with any of the solutes under inves- 
tigation in this work increases the equilibrium volume 
(Table II) . Unconstrained uniaxial tensile loading is also 
accompanied by an increase of the volume per atom, as 
depicted in Fig. 7 for Fe and for the Fe-based alloys with 
the highest considered solute concentrations. In the re- 
gion of the maximum strains, both Fe and the Fe alloys 
have rather similar volumes as opposed to the situation 
at equilibrium. The volume change of Fe due to the 



is at most 1.5% (for one and the 
^ is quite similar to e,n for Fc and 



deformation is thus distinctly larger than those of the 
Fe-based alloys. According to our numerical data (also 
Fig. 6), the larger volume increase of Fe is paralleled by a 
larger magnetic moment increase along the uniaxial de- 
formation from equilibrium to e^n, despite the fact that 
Fe contributes only by a factor of 1 — a; to the magnetic 
moment of the alloys. 

The comparatively large increase of the magnetic mo- 
ment in Fe along the uniaxial deformation path with un- 
constrained volume and the pronounced difference be- 
tween cr^ and fJni (5 GPa) may, at least in parts, be 
related. We recall, that both the increase of the mag- 
netic moment and the difference a^ — a^ of all Fe-alloys 
are smaller than the values of Fe. To substantiate this 
correlation, we recalculated the ITSs for the previously 
determined unconstrained strain paths but with fixed 
magnetic moments, i.e., the magnetic moments were not 
allowed to relax to self-consistency. We employed the re- 
spective ground state magnetic moments of Fe and the 
Fe-V alloys. 

Fixing the magnetic moment increases the resulting 
computed strengths, but the increase is strongly dimin- 
ished with higher solute concentration, e.g., approxi- 
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FIG. 7. (Color online) The volume per atom of Fe and of 
Fe-based alloys with 10 % solute concentration (5 % solute 
concentration in the case of Mn) as a function of strain for the 
unconstrained deformation path. The open symbols specify 
the volumes at the maximum strain. 



mately 4.0 GPa for Fe but only 0.2 GPa for Fe-lOV. 
The previous numbers may be obtained by comparing (Tm 
from Fig. 1 with the ITSs obtained with the fixed-spin 
moment constraint depicted in Fig. 4 (open pentagons). 
The latter strength values are plotted in comparison to 
the previously determined ITSs computed without a con- 
straint on the magnetic moments but for fixed volume 
(ct^). The alloying trend that both curves for Fe-V fol- 
low in Fig. 4 is alike, and the resulting stresses are es- 
sentially much larger than cr,n of Fe (4.0 GPa for the 
fixed-spin moment value and 5.0 GPa for a^J but only 
slightly to moderately larger than a^ of Fe-lOV (0.2 GPa 
for the fixed-spin moment value and 1.6 GPa for cj^). 
Thus, fixing the volume but allowing for a relaxation of 
the magnetic moments produces essentially the same al- 
loying effect than fixing the magnetic moment but taking 
into account structural relaxations. 

However, as we pointed out above, the magnetic mo- 
ment is almost unchanged along the constant volume de- 
formation owing to the weak influence of the c/a ratio. 
Thus, we conclude, that the main contribution to the 
difference of a^^ and cr,n originates from a magnetic con- 
tribution to the total energies which is associated with 
the relaxation of the spin magnetic moment to its self- 
consistent value. This magnetic contribution is accord- 
ingly the largest for Fe, which can be partly ascribed to 
the comparatively small volume at equilibrium. Both the 
reduced volume change and the involved reduced mag- 
netic contribution to the stress bring forth the conver- 
gence of a^i to (Tin for the investigated Fe-alloys with 
10% solute concentration. 

Returning in conclusion to Eq. (3), the correlation be- 
tween cr,n and cr^^ fails due to magnetism in Fe and 
the investigated Fe alloys. That is, the alloying trend 
on the fcc-bcc SEDs at fixed volume can not capture the 



alloying trend on the energy-strain landscape and hence 
the ITS of the uniaxial deformation path as opposed to 
nonmagnetic V-Ti-Cr alloys. ^^ 



V. CONCLUSIONS 

The ideal tensile strengths of bcc ferromagnetic Fe and 
Fei-xM^ (M=Cr, Ni, Al, Co, Mn, and V) random alloys 
in the [001] direction have been investigated using the all- 
electron exact muffin-tin orbitals method in combination 
with the coherent-potential approximation. The present 
calculated ideal strength value of Fe and the branching 
point from the primary bet deformation path to the sec- 
ondary orthorhombic deformation path of Fe agree well 
with previously published results, which thus confirms 
that our methodology has the accuracy needed for such 
kind of calculations. For the Fe-based alloys, we found 
that the ideal strength increases with increasing concen- 
tration of Cr, Co, and V and decreases with Ni and Al 
addition into pure Fe. Mn shows a non-monotonous al- 
loying behavior. Unlike the nonmagnetic bcc V-based 
alloys investigated in our previous paper, ^■^ constant vol- 
ume fcc-bcc structural energy differences can not entirely 
capture the alloying effect on the ideal tensile strength 
for the FM bcc Fe-based alloys. By calculating auxiliary 
ideal strengths assuming a constant volume deformation 
(cr^) and fixed-spin moment calculations along the (kept 
unaltered) tetragonal deformation paths, we showed that 
mainly the interplay between the volume increase and 
magnetic moment increase along the tetragonal defor- 
mation paths lead to a failure of the correlation between 
am and af^^. 

The present results offer a consistent starting point 
for further theoretical modeling of the micro-mechanical 
properties of Fe-based alloys. Based on these achieve- 
ments, we conclude that the EMTO-CPA approach pro- 
vides an efficient and accurate theoretical tool to de- 
sign the mechanical strength of ferromagnetic bcc ran- 
dom solid solutions Nevertheless, in such applications one 
should always monitor the basic muffin-tin and single-site 
GPA errors and make sure that they remain at acceptable 
level as a function of the lattice distortion and chemical 
composition. 
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